Multi-photon spectroscopy of a hybrid quantum system 
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We report on experimental multi-photon spectroscopy of a hybrid quantum system consisting of a 
superconducting phase qubit coherently coupled to an intrinsic two- level defect. We directly probe 
hybridized states of the combined qubit-defect system in the strongly interacting regime, where 
both the qubit-defect coupling and the driving cannot be considered as weak perturbations. This 
regime is described by a theoretical model which incorporates anharmonic corrections, multi-photon 
processes and decoherence. We present a detailed comparison between experiment and theory and 
find excellent agreement over a wide range of parameters. 



I. Introduction 

One of the fundamental tenets of quantum mechanics 
is that the rules that govern a system's behavior do not 
depend on its physical size. Rather, it is the number 
of degrees of freedom a system possesses and how well 
the system is isolated from its environment, which dic- 
tates how it behaves and whether it is considered a quan- 
tum or classical system. This is most clearly seen in the 
plethora of different systems in which quantum mechan- 
ical phase coherence can be seen and controlled. These 
include atomic (sub nm) systems such as trapped ions^ 
or optical color-centers^, mesoscopic systems (nm-/im) 
such as superconducting devices^ or macroscopic objects 
(mm) such as Bose-Einstein condensates^. 

In this paper, we report on experiments with a hy- 
brid quantum systenP consisting of two parts of different 
length-scale and physical nature, namely a superconduct- 
ing phase qubit (~ /im) and an intrinsic two- level defect 
(~nm). The systems are coupled via an effective dipole- 
dipole interaction in a regime where both the interaction 
between its components and the driving field are strong. 
We observe effects which are not present in the usual sim- 
plified models^, including multi-photon transitions and 
weak violation of selection rules due to higher order cor- 
rections to the anharmonic potential. 

When studying superconducting qubits spectroscopi- 
cally, one often finds clear signatures of level anticrossings 
at certain frequencies^^i^. These anticrossings indicate 
intrinsic, microscopic two-level systems (TLS) being cou- 
pled coherently to the qubit circuit. In general, ensembles 
of two-level microscopic defects are believed to be respon- 
sible for loss in a wide variety of systems, including phase- 
and flux-based superconducting qubits and even nanome- 
chanical oscillators^^ ^^, as well as more general effects in 
amorphous and spin-glass systems^^. In superconduct- 
ing systems, the exact microscopic nature of these defect 
two-level systems remains unknown. It has been shown 
that they are coherent two- level systems, and due to their 
relatively long coherence times can potentially be used 
to store and retrieve quantum information^. In gen- 



eral, two-level defects are considered detrimental to the 
qubits operation, since they introduce additional chan- 
nels of decoherence^^ ^^ and therefore, acquiring a more 
in-depth understanding of their nature is essential. It has 
been suggested that these TLSs are formed by charged 
microscopic defects, located inside the insulating tunnel 
barrier of the Josephson junction^ ^^ and interacting with 
its electric field. In this case, the coupling to the qubit 
circuit can be written as ex d/ • E, where df is the dipole 
moment of the TLS and E is the electric field across 
the Josephson junction (JJ). Another suggested mech- 
anism for qubit-defect coupling is that the state of the 
TLS modulates the transparency of the junction^ ^^ and 
therefore its critical current Ic- In this case two- level 
defects could be form ed by Andreev bound states^ or 
Kondo-impurities^^^. Since the form of the potential in 
the qubit circuit depends on the strength of the critical 
current, a modulation of Ic will always result in a direct 
change of the level splitting in the circuit. Since we do 
not observe signatures of longitudinal coupling in our ex- 
periments (as has also been shown in flux qubits^^), we 
rule out this coupling mechanism and take the TLS to 
be a charged defect. 



II. Experiment 

We measured a superconducting phase qubit, consist- 
ing of a single Josephson junction embedded in a super- 
conducting loop (Fig.[T]). The experiments reported here 
are performed with the qubit circuit biased such that the 
frequency of its lowest lying transition is close to the res- 
onance frequency of a particular TLS. The studied TLS 
defect has a level splitting of Cf/h = 7.845 GHz and is 
coupled to the qubit with a coupling strength v±/h = 21 
MHz. 

The J J is shunted with a capacitor Cs of 800 fF, which 
is needed to tune the circuit to the desired frequency 
range. Its critical current is Ic = 1.1 /iA and the induc- 
tance of the loop is L = 720 pH. During the experiment, 
the chip is maintained at a temperature of approximately 
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FIG. 1: Sketch of the experimental setup and circuit of the 
phase qubit with the TLS residing inside the qubits Josephson 
junction. 



40 mK. 

By applying an external magnetic flux through the 
loop, the shape of the potential of the quantum circuit 
can be tuned. In the parameter regime of interest, the po- 
tential consists of two wells, one of which is very shallow. 
Operation of the circuit takes place in this anharmonic 
shallow well, containing five to ten energy levels^. The 
two lowest energy states in this well form a macroscopic 
two-level system (our qubit). The transition frequencies 
are in the microwave range and the state of the system 
can be manipulated via capacitively coupled microwaves. 
In order to readout the state of the qubit, the potential 
barrier between the two wells is lowered to allow the ex- 
cited states of the shallow well to tunnel into the deeper 
well. States in the two wells differ by a large amount 
of persistent current through the loop, which can be de- 
tected by measuring the switching current of a measure- 
ment SQUID located on the same chip. The resulting 
escape probability is normalized to the measurement con- 
trast defined by the qubit-SQUID coupling. A measure- 
ment thus distinguishes between the ground state and all 
other, excited states of the qubit circuit and yields the 
overall excitation probability for the qubit or the proba- 
bility not to be in the ground state. 

In the following we present data on microwave spec- 
troscopy of the circuit. The system is driven for a suffi- 
ciently long time, ~ 2 /is, so that all transient phenomena 
have decayed. The state of the qubit after this pulse is 
measured as a function of microwave frequency and fiux 
bias to map out the system transitions. 



III. Theory 

When considering a phase qubit which is weakly driven 
and/or weakly coupled to a TLS, it is common to model 
the qubit circuit as a two-level system and to use the 
rotating- wave- approximation. For our experiment, nei- 
ther of the above approximations are valid and one must 



use a formalism which takes into account an arbitrary 
number of qubit states as well as multi-photon processes. 
This allows us to describe the experiment over a wide pa- 
rameter range without knowing which effects are impor- 
tant a priori. Several techniques exist for treating a driv- 
ing field when the perturbation is no longer "weak''^^^^^. 
We employ a numerical expansion in the Floquet ba- 
sis^^ 27 g^g |-]^^g allows us to include arbitrary multi-photon 
processes together with decoherence. The inclusion of 
higher lying states in the qubit can be achieved via di- 
rect diagonalization of the exact qubit potential, or by 
using an anharmonic approximatioiPEU^ 

The qubit is well described by a quantum anharmonic 
oscillator model, with the anharmonicity determined by 
the Josephson energy Ej^. In the parameter regime 
explored in the experiment, it is sufficient to describe 
the qubit by using only the lowest three energy levels. 
The TLS is modeled as a two-level system with its en- 
ergy splitting Cf obtained from the spectroscopic data. 
We write the coupling term between qubit and TLS as 
v_Lp-rx^ where the operator p is proportional to the charge 
and therefore the electric field across the qubit circuit's 
Josephson junction. The operator Tx is the Pauli matrix 
in x-direction acting on the two-level defect, inducing 
transitions of the state of the TLS. The Hamiltonian is 
given by 

^0 = eoi |1) (l| + (eoi+ei2) |2) (2| + -6/f, + -^xp-r, . (1) 

The operator p can be written in the harmonic oscillator 
basis as oc i{d^ — a) ^ where a, a^ are the usual annihila- 
tion and creation operators. In the exact qubit eigenbasis 
|n), used to write down Eq. (fTl), this operator will acquire 
corrections of the form |n) (n + 2|. To obtain the energy 
splittings eoi and ei2 in this basis, we solve numerically 
using the exact potential of the phase qubit circuit. The 
state of the TLS is described by its ground state \g) and 
excited state |e). The energies are defined relative to the 
energy of the ground state of the coupled qubit-TLS sys- 
tem |0^). The system is coupled to an external microwave 
drive at frequency ujd with the coupling Hamiltonian 



^i = -Vqp cos Udt, 



(2) 



where the coupling constant Vq is proportional to the 
microwave amplitude and we define the generalised n- 
photon Rabi frequency ftq oc V^ on resonance. 

Previously, we have shown that close to resonance, the 
TLS is driven by a quasi-direct driving term mediated by 
a second order process involving the qubit circuit's third 
leveP^. As the model here already includes the third level 
by construction, this effect is incorporated automatically. 

To describe the effects of decoherence in the system, 
we use a simple Lindblad modeP^ ^^, with relaxation and 
dephasing rates taken from independent measurements 
of both qubit and defect. In numerical modeling of the 
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system, we use the values of T^^ = 120 ns and T{ 
715 ns for the relaxation rates of the qubit and defect. 



respectively. Dephasing rates , determined from Ramsey 
fringes, are found to be T2 = 90 ns for the qubit and 

T2 = 110 ns for the TLS. The dynamical equations for 
the reduced density matrix of the system (consisting of 
qubit and TLS) are given by 
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where the sum is over all possible channels of decoher- 
ence with the respective rates Vj. The Lj are the op- 
erators corresponding to each decoherence channel, e.g., 
pure dephasing of the qubit is described by the operator 
a) a while relaxation of the qubit is described by the ani- 
hilation operator a. This formalism in principle enables 
us to include an arbitrary number of qubit levels without 
changing the structure of the theory. 



IV. Results and Discussion 

A. Two-photon spectroscopy of 4-level hybrid 
quantum system 



We initially present results from driving the qubit with 
a relatively low microwave power (Fig. E|, corresponding 

to a qubit Rabi- frequency oi Viq /h ^ 4 MHz. In this 
regime, it is sufficient to approximate the qubit circuit 
by a two-level system and we can write the Hamiltonian 
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with Pauli matrices for the qubit a and for the defect r. 
A sketch of the level structure including the relevant tran- 
sitions is shown in the inset of Fig. [2] Here we define the 
hybridized eigenstates |1 — ) = cos^|Oe) — sin^|l^) and 
|1-F) = sin6> |Oe)+cos^ |1^), where tan 2^ = 'y±/(eoi — e/) 
varies with flux bias. The coupling term cFyTx is formally 
equivalent to the previously described dipole-dipole cou- 
pling ~ E ■ df. The experimental data in Fig.pfa) shows 
a characteristic level anticrossing associated with the 
qubit-TLS resonance, eoi ~ e/, with coupling strength 
v^/h = 21 MHz. 

An additional spectroscopic line can be seen in the 
middle of the qubit-defect anticrossing. We identify this 
line as a two-photon transition from the ground-state |0^) 
to the excited state |le), shown by red arrows in the inset 
of Fig. [2] Similar two-photon absorption mediated by 
dipole-dipole coupling is observed in experiments with 
pairs of dye molecules at high excitation powei^^. 

In Fig. [2{b) we show the result of a theoretical treat- 
ment of the system, calculating the time dependence of 
the system's density matrix including decoherence. Sys- 
tem parameters are taken from the measurements as de- 
scribed above. As can be seen, the numerics reproduce 
the experimental data with high accuracy. 
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FIG. 2: (Color online) Experimental (a) and theoretical (b) 
spectroscopic scans as a function of flux bias and drive fre- 
quency, showing both the anticrossing due to single-photon 
transitions and the resonance line associated with the weaker 
two-photon transition. The colorscale indicates the probabil- 
ity to find the qubit not in the ground state. We use the 
same colorscale for both plots. Inset: Energy level structure 
of the qubit/defect system including the hybridised states 
|1±), which are split by the coupling v±. As well as the 
usual single-photon transitions (between lO^f) and |1±)), a 
two-photon transition is allowed between the states |le) and 



The power dependence of the two-photon feature is 
shown in Fig. [3J to confirm that the observed middle line 
results from a two-photon transition. In this figure the 
amplitudes of the two-photon (|0^) ^ |le)) and the lower 
single-photon (|0^) ^ |1— )) transitions are plotted as a 
function of driving power. The data is fitted with the 
conventional saturation curve^^ 

[i + (!^(")/r)2]' ^^^ 

where Q^^'' /T is proportional to the ratio of (n-photon) 
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FIG. 3: (Color online) Detailed power dependence of both the 
single (top, solid blue) and two-photon (lower, dashed red) 
transitions close to resonance with the defect (as shown in 
Fig. [2|. The amplitudes and associated uncertainties are ob- 
tained by fitting the experimental data with Lorentzian func- 
tions. The scaling of the saturation curves is as expected for 
both single- and two-photon transitions, as indicated by the 
fitted lines (see text). 



Rabi frequency to loss rate and n = 1, 2 is the order of the 
photon transition, confirming the single and two-photon 
nature of these two transitions, respectively. 

The position of the two-photon line in the middle of the 
anticrossing gives an indication of the strength of a possi- 
ble longitudinal coupling ~ a^f^ between qubit and TLS. 
Through the presence of the two-photon line exactly in 
the middle of the anticrossing, we can conclude (i) the 
TLS is a two-level or at least strongly anharmonic sys- 
tem and (ii) longitudinal coupling between qubit circuit 
and TLS is negligible. Very similar features have been 
observed earlier in a related system, namely a supercon- 
ducting flux qubit ^^, which suggests that these strongly 
coupled TLS have the same origin in both flux and phase 
qubit s, even though the degrees of freedom being manip- 
ulated are different. Additional evidence against the TLS 
being an harmonic object (as would result from an en- 
semble of two- level systems forming an effective TLS^^, 
for example) is provided by experiments trying to pump 
two excitations resonantly into the TLS^^, similar to the 
protocol used in Ref. [16i 



B. Two-photon spectroscopy at strong excitation 



In addition to the structure around the defect-qubit 
anticrossing, at higher power {Qq ^ /h ^ 7 MHz) we find 
additional features at lower flux bias, shown in Fig. |4] 
We observe a weaker line running parallel to the single- 
photon qubit transition, showing an additional anticross- 
ing with the |0^) ^ |le) line, indicated by the lefthand 
dashed box on the figure. From their power dependence, 
both these lines can be identified as two-photon transi- 
tions. 
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FIG. 4: (Color online) At higher microwave powers and over 
a wider region of parameter space, processes involving higher 
lying states of the qubit become important. Experimental (a) 
and theoretical (b) spectroscopic scan as a function of flux 
bias, showing both the transitions to the hybridised states 
|1±) (as per Fig.lSl and the states |2zb). The colorscale is the 
same for both graphs. Inset: Level structure including higher 
lying states of the qubit including further hybridised states 
which can be excited via a two-photon process, as indicated. 



To explain the spectral features of Fig. [4ja), it is now 
essential to include the third level of the qubit circuit in 
the treatment. A sketch of the level structure around 
this second anticrossing is again shown in the inset of 
Fig. |4] Possible transitions are indicated by arrows. Due 
to the anharmonicity in the qubit circuit, the transition 
between the levels |1) <-> |2) is detuned from the qubit 
transition |0) ^ |1) by the amount A = ei2 — eoi < 0. 
At the chosen range of flux bias, our circuit has an an- 
harmonicity |A| ~ 100 MHz. We therefore identify this 
weaker, parallel line as the two-photon |0^) ^ \2g) transi- 
tion, detuned from the first transition line by the amount 
I A 1/2 ~ 50 MHz due to its two-photon nature. The an- 
ticrossing in the left dashed box is due to the coupling 
between qubit and TLS ~ P ■ Tx^ which hybridizes the 
states |le) and \2g). The coupling strength is increased 



by a factor \f2 due to the properties of the momentum 
operator and the magnitude of the anticrossing is given 
by V^'U^/2 ^15 MHz (where the additional factor of two 
results from the two-photon process). Fig.[4Fb) shows the 
result of a numerical calculation. The theory is in very 
good agreement with the experimental data over all pa- 
rameter regimes. 

An additional feature can be identified in the spec- 
troscopic scan of Fig. [41 a) at a flux bias of 1.175 and 
microwave frequency resonant with the TLF ujd ^ ^f ^ 
^602- Here, the two-photon |0^) ^ \2g) line crosses the 
one photon |0^) ^ |0e) transition. This is an area of high 
symmetry in the spectrum as is illustrated in Fig. [5J re- 
sulting in strongly enhanced absorption. The asymmetric 
nature of this effect is well reproduced by the numerics 
(Fig.gb)). 




FIG. 5: (Color online) Level structure of the coupled system 
of qubit and TLS for ujd ^ ef ^ |eo2, illustrating the high 
degree of symmetry. Possible microwave induced transitions 
at this driving frequency are indicated by arrows. 



C. Single photon spectroscopy of upper hybridised 
states 
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FIG. 6: (Color online) Using higher frequency excitations, 
the hybridised states |2di) can be directly excited via a one- 
photon process. Experimental (a) and theoretical (b) spec- 
trum showing the clear anti-crossing due to the lO^^) ^ |2zb) 
transition. The colorscale is the same for both graphs. (In- 
set) The corresponding energy level diagram shows this direct 
transition which is typically ignored within the usual simpli- 
fied models of such a system. 



As a further confirmation that the left anticrossing of 
Fig. 4 indeed represents the coupling between the states 
|le) and \2g), we probed these higher lying states directly 
with a single-photon transition. Such a process is equiv- 
alent to the cooperative absorption of single photons by 
pairs of rare-earth ions in crystals^^. The experimental 
results are presented in Fig. 6ja) where one can clearly see 
the anticrossing at doubled resonance frequency ~ 15.7 
GHz with a magnitude now of \f2v\_ ~ 30 MHz. When 
approximating the shallow well using the third order an- 
harmonic oscillator potential^ the transition |0) ^ |2) 
is forbidden. However, when treating the anharmonicity 
exactly, there is a finite matrix element in the momen- 
tum operator p corresponding to this transition. Fig.JGJ^b) 
shows the results of a theoretical calculation and the inset 
gives a sketch of the level structure including the transi- 
tions induced by the microwave drive. 



V. Conclusion 

In conclusion, we have presented detailed data on ex- 
periments which probe the hybridization between a su- 
perconducting phase qubit and a two-level defect in the 
tunnel barrier of its Josephson junction. These ele- 
ments are all strongly interacting, forming a regime anal- 
ogous to strong coupling in circuit QED, with the role 
of the photon mode played by the TLS. The resulting 
hybrid quantum system is in turn probed with multi- 
photon spectroscopy. Understanding the experimental 
results requires a theoretical model which must go be- 
yond the standard weak driving, harmonic well and two- 
level approximations. Using this tight link between preci- 
sion spectroscopy and systematic theoretical models, we 



are able to map out the quantum mechanical nature of 
this hybrid system, combining elements from two length 
scales differing by four orders of magnitude. 
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